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1 Introduction 



As the LHC begins operation, the number one topic of the day is QCD. To a 
first approximation, everything that happens at the LHC is determined by QCD. 
Events from new physics at the LHC will appear, if we are lucky, at fractions of a 
part per billion of the proton-proton total cross section. They will appear at rates 
ten thousand times smaller than Drell-Yan processes, and a hundred times smaller 
than top quark pair production. This requires that we understand those Standard 
Model processes at a level that includes their dressing with QCD radiation. 

Since the 1980's, there has been considerable work on methods for QCD com- 
putation. The primary goal of this work has been to make seemingly impossible 
calculations doable. In the early 1980's, it was beyond the state of the art to com- 
pute to lowest order the rates at hadron colliders for production with three jets. 
Recently, this process has been computed at the 1-loop level, a computation that, 
if done by textbook methods, would have required thousands of Feynman diagrams, 
each one an integral over millions of terms [HE]- 

A byproduct of this improved understanding of QCD computation is that calcu- 
lations of reasonable difficulty by textbook methods become trivial when approached 
with the new methods. In the LHC era, every graduate student ought to be able to 
calculate the QCD amplitudes for multijet processes. In this lecture, I will give you 
some tools to do that. 

I will first review a set of ideas developed in the 1980's. Sections 2, 3, and 4 will 
develop, in turn, the ideas of spinor products, color ordering, and MHV amplitudes. 
In Section 5, I will illustrate how these methods make the somewhat strenuous cal- 
culation of the most basic quark and gluon cross sections a triviality. In Section 6, 
I will review a relatively new wrinkle in this technology, the Britto-Cachazo-Feng- 
Witten recursion formula. In Section 7, I will apply that formula to compute the 
Drell-Yan cross sections with one and two jets. I hope that this last example will give 
a persuasive illustration of the power of these methods. 

There are many references for those who would like to learn more. The first three 
topics are reviewed very beautifully in a 1991 Physics Reports article by Mangano and 
Parke |3] and in the 1995 TASI lectures by Dixon [4J. The newer set of calculational 
tools are described in recent review articles by Bern, Dixon, and Kosower [5] and 
Berger and Forde [B]. These last articles are especially concerned with technology 
that extends the ideas presented here to the computation of loop amplitudes. 

Loop calculations, though, are beyond the scope of these lectures. My only goal 
here is to do easy calculations. I hope that these lectures will help you extend the 
domain of QCD calculations to which you will apply this term. 
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In parallel with learning these methods, it will be useful for you to explore some 
of the available computer programs that automatically generate matrix elements 
for multi-particle processes and also carry out Monte Carlo integration over phase 
space, a topic that is not discussed in these lectures. A particularly accessible and 
powerful code is MADGRAPH/MADEVENT [HH]. An important advanced code is 
SHERPA piTO] with the matrix element generator COMIX [11]. It is always good, 
though, to know what is inside the black box and to have tools for extending what is 
available there. For that, I hope these lectures will provide useful background. 



2 Spinor products 

Perturbative QCD is primarily concerned with the interactions of gluons and 
quarks at momentum scales for which the masses of these particles can be ignored. It 
is well-known, and discussed in the textbooks, that massless particles can be labeled 
by their helicity. For massless states, the helicity is a well-defined, Lorentz-invariant 
quantity. The basic goal of these lectures will be to compute tree amplitudes for 
massless quark and gluon states of definite helicity. 

In the 1980's, Berends and Wu spearheaded an effort to compute amplitudes for 
massless particles by exploiting the property that their lightlike momentum vectors 
can be decomposed into spinors [T^. Using this idea, scattering amplitudes can be 
written in terms of Lorentz-invariant contractions of spinors, spinor products. 



2.1 Massless fermions 



Consider a massless fermion of mommentum p. The spinors for this fermion satisfy 
the Dirac equation 

^ U{p) = (1) 

There are two solutions to this equation, the spinors for right- and left-handed 
fermions. In a basis where the Dirac matrices take the form 

where a'^ = (1, a), = (1, —(f), these spinors take the form 

where the entries are 2-component spinors satisfying 

p. cr ur = , p-auL = 0. (4) 
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Each equation has one unique solution. The spinor UFt{p) can be related to the spinor 

uiip) by 

ur(p) = icr'^ulip) , (5) 

since this transformation turns a solution of one of the equations Q into a solution 
to the other. This equation also gives a phase convention for ur{p) that I will use 
throughout these lectures. 

To describe antiparticles, we also need solutions V{p) that describe their creation 
and destruction. However, for massless particles, V{p) satisfies the same equation as 
U{p). We can then use the same solutions for these quantities, with ur and ul related 
as above. The spinor Vr{p) is used for the creation of a left-handed antifermion; the 
spinor Vl{p) is used for the creation of a right-handed antifermion. 

In the discussion to follow, it will be simplest to treat all particles as final states 
of the amplitudes we are considering. The outgoing left- and right-handed fermions 
will be represented by the spinors Ul and Ur and outgoing left- and right-handed 
antifermions will be represented by the spinors Ur and Ul. I will recover initial-state 
particles by crossing, that is, by starting from the amplitude with the antiparticle in 
the final state and evaluating that amplitude for a 4-momentum with negative time 
component. 

I will represent these spinors compactly as 

Ul{p) = {p , Ur{p) = [p , ULip) = p] , Ur{p) = p) , (6) 

The Lorentz-invariant spinor products can then be standardized as 

UL{p)UR{q) = {pq) , UR{p)UL{q) = [pq] . (7) 

The quantities on the right-hands sides are called simply angle brackets and square 
brackets. The spinors are related to their lightlike 4-vectors by the identities 

P)[p = URippRip) =^(^^) , P]{P = Ul{p)Ul{p) =A^Af-) (8) 
From these formulae, we can derive some basic properties of the brackets. First, 

{pq) = [qp\* . (9) 
Next, multiplying two matrices from ([s]) and taking the trace, 

{pq)[qp\ = ii[4 A^A^)] = '^p-q (10) 

so that 

\{pq)\'=\[qp]\' = 2p-q. (11) 
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Finally, using (|5|, 

(pq) = u{{p)uR{q) = ul^{p){ia^)abul^{q) ■ (12) 
Then, by the antisymmetry of cr^, 

{pq) = -{qp) [pq] = -[qp] . (13) 

We now see that the brackets are square roots of the corresponding Lorentz vector 
products, and that they are antisymmetric in their two arguments. 

Some further identities are useful in discussing vector currents built from spinors. 
First, 

ULiph'ULiq) = u[{p)a>'uLiq) , Uniph'URiq) = ui{p)a>'uRiq) . (14) 
The relationship between ui{p) and uji{p) in ^ implies that we can rearrange 

u[{pp^UL{q) = u{{p)a^{-{zay)uL{q) 

= ul{p){-zcr')cr^^{za')u,{q) 

= ul{p)a>^^uUq) 

= u^R{q)a^UR{p) , (15) 

From this, it follows that 

{pi^q] = [qi^p) ■ (16) 

The Fierz identity, the identity of sigma matrices 

{(r^)abicf^),d = 2{ia'),,{ia')bd , (17) 
allows the simplification of contractions of spinor expressions. In a matter similar to 



(15), the relation (17) can be used to show 



(P7^g] {k7,e] = 2{pk) [iq] , (p7^g] [kl,i) = 2{pi) [kq] . (18) 

Finally, spinor products obey the Schouten identity 

{^J){k£) + {^k){£J) + {^i){Jk) = 

[tj][ki] + [ik][ij] + [ti][jk] = . (19) 

To prove this identity, note that the expressions on the left are totally antisymmetric 
in j,k,£. But, antisymmetrizing three 2-component objects gives zero. 

The identities presented in this section will allow us to reduce complex spinor 
expressions to functions of angle brackets and square brackets, which are the square 



7 



Figure 1: Feynman diagram for e^^e^ — )• /U^^/uJ. 



roots of the Lorentz products of lightlike vectors. In the following, I will denote these 
Lorentz products using the notation 

Sij = 2pi ■ pj = (pi +Pjy . (20) 

To evaluate simple expressions written in terms of spinor products, we can often 
think of the whole expression as the square root of an expression in terms of the Sij. 
However, for more complex expressions, it is usually easiest to directly evaluate the 
spinor products from their component momenta. I give formulae for evaluating spinor 
products in Appendix A. 



2.2 e+e 

At this point, we are already able to perform some interesting computations. 
Consider, for example, the tree-level amplitude for e^e^ — )■ f^Lfi~ji in QED. This 
is given by the diagram shown in Fig. [T} Label the momenta as in the diagram, 
considering all momenta as outgoing. Then the amplitude is 

^M = i-tef ^ f7z.(3)7^f/L(4)f7z.(2)7^f/i(l) 
= -(37'^4](27,l] 

= ^(32) [14], (21) 



where I have used (18) in the last step. Now, (32) and [14] are both square roots of 

S23 = {k2 + hf = {h + hf (22) 

which is just the Mandelstam invariant u. In the e^e~ center of mass frame, u = 
-2^2(1 + cos^), and = s = AE^. Then 

\tM\^ = e\l + cosey . (23) 
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This is a familiar result, and we have obtained it with surprising ease. 



Actually, we can simplify the result further. In the denominator of (21), q = 
2ki ■ k2 = (12) [21]. Multiply the numerator and denominator of (21) by (32), to give 

_ ,,,.(32)1141(32) 

(12) [21] (32) ^ ^ 

Then, in the denominator 

[21] (32) = [12] (23) = [1 ? 3) = [1(- /- ^- ^)3) (25) 

Using the Dirac equation, = 0, ^3) = 0. This leaves (—[14] (43)). The last square 
brackets cancel and we find 

tM=2te'-^^f^ . (26) 

(12) (34) ^ ' 

The entire expression can be written in terms of angle brackets with no square brack- 
ets. Similarly, if we had multiplied instead by 

M (27) 
[14] ^''^ 

a similar set of manipulations would have given 



2 [14] 
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= MSI • 



with square brackets only. 



2.3 Massless photons 

This simplification of amplitudes with fermions extends to amplitudes with mass- 
less vector bosons. I will show that the polarization vectors for final-state massless 
vector bosons of definite helicity can be represented as p^l4pi5] 

Here k is the momentum of the vector boson and r is some other fixed lightlike 4- 
vector, called the reference vector. The only requirement on r is that it cannot be 
coUinear with k. 



9 



Here is the argument: First, note that (29) satisfies the basic properties 

The second condition follows from ^k] =0. Also, 

■ [^im* = eUk) ■ e^k) = , 



because, by (18), this expression is proportional to (rr) = 0. Finally, 

I * /r N|2 ^ 1 jrYk] {k-j^r] ^ 2 {rk) [rk] ^ _ 
2 {rk)[kr] 2{rk)[kr] 

-1 , e^(A:) ■ [emr = • 



Then 



i4(^)r = ki(^)r 

as required for polarization vectors. 



(30) 
(31) 

(32) 
(33) 



Next, evaluate the formulae for the e*^^^ for a particular choice of the reference 
vector r. For k = {k, 0, 0, k), let r = (r, 0, 0, — r). The associated spinors are 



UL{k) = V2ki^j UR{k) = V2k 



Uiir) = v2r 



(34) 



Then 



{rYk] = ulirpf'uLik) = VAk^ (0, 1, -i, 0)^ 

(rk) = u\{r)uL{k) = —a/ 4kr (35) 

For this choice of r, e*j^{k) is manifestly the right-handed polarization vector. 



4'^(fc) = -^(0,1,^,0) 



(36) 



Finally, analyze how (29) changes when we change the reference vector r to a 
different lightlike vector s. 



t^{k-r)-ej!^{k-s) 



1 UrYk] {sYk]\ 

l.^^i-irrk]{ks)Hsrk]{kr)) 

^^(-H'^^.) + (.7^^r)) 

1 1 
^/2{rk){sk) 
1 1 

^{rk){sk) 
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■ 2A;^(sr) . 



(37) 



Figure 2: Ward identity obeyed by a gauge-invariant sum of diagrams with all external 
particles on shell. 




The last line follows from the anticommutator of Dirac matrices. The final result of 
this calculation is that 



e 



^^(A;;r)-e*/(fc;s) = /(r,s)P. (3^ 



where /(r, s) is a function of the two reference vectors. This expression, when dotted 
into an on-shell photon or gluon amplitude, will give zero by the Ward identity, as 
shown in Fig. [2j Thus — as long as we are computing a gauge-invariant set of Feynman 
diagrams — we can use any convenient reference vector s and obtain the same answer 



as we would with the particular reference vector r used in (36). 



This completes the justification of the representation (29) of massless photon or 
gluon polarization vectors. From here on, because I will in any case consider all 
momenta as outgoing, I will drop the explicit *'s on the polarization vectors. 



2.4 e'^e — )■ 77 



We can illustrate the application of these polarization vectors by computing the 
amplitudes for e^e^ — > 77. Label the momenta as in Fig. [3| taking all momenta as 
outgoing. Then the value of the amplitude for this process is 

= (-26)2(2 (7 ■ e(4)^^^^7 ■ e(3) + 7 ■ e(3)^^^^7 ■ e(4)| 1] . (39) 
I S24 S23 I 
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We will choose explicit polarization vectors e(3), e(4) to evaluate this amplitude for 
each choice of photon helicities. This formula also introduces some additional stream- 
lining of the notation: I write (2 + 4), instead of {^2+ ^4) or even (^+ and I use 
(20) to express the denominators. 

There are four possible choices for the photon polarizations. However the cases 
7k7l and 7l7_r are related by interchange of the momenta 3 and 4. The cases 7k7r 
and 7l7l are related by parity, which interchanges states with R and L polarization. 
Further, it is easy to see that the amplitudes for •jR'yR and 7l7l are actually zero. 
For the case of 'Jr'Jr, choose r = 2 for both polarization vectors, 

^ ' V2 (23) ' ' ' 72 (24) ' ' 



When these choices are used in (39), we find, with the use of the Fierz identity (18) 

(27-6(4) ~ 2(22) [4 = , (41) 

which vanishes because (22) = 0. A similar cancellation occurs with e(3). So the 
entire matrix element vanishes. The amplitude for the case 7l7l must then also 
vanish by parity; alternatively, we can find the same cancellation for that case by 
using r = 1 in both polarization vectors. 

To compute the amplitude for the case ■JrJl, choose 

6^3) = e-(4) = (42) 

V2 (23) ' V2 [14] • ^ ' 

Then the second diagram in Fig. |3] vanishes by the logic of the previous paragraph. 
Using the Fierz identity, the first diagram gives 

—ie^ 2 ■ 2 
'-^ = 1^ (-2)(23)|141 <^^>I^P + ^>^>P'I 

(24) [14] (42) [31] 



si3(23)[14] 
(13)[311(23)|14] <^^>'^^i<^^>P^' 

2^e'^iM)l (43) 
{23)(31) ^ ' 



In terms of the Mandelstam variables, S23 = u, S13 = S24 = t, so 



|,M|2 = 4e^l =4e^L_^ , (44) 
' ' u 1 + cos 6^ ^ ^ 
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Figure 4: Feynman diagrams for qlQr — ^ 99- 



which is the well-known correct answer. In this case, the calculation gave directly 
the simplified form of the expression that involves only angle brackets and no square 
brackets. 

The methods of this section can be generalized to massive external particles. For 
the treatment of W and Z bosons, see Section 7. An efficient formalism for the 
treatment of massive fermions is presented in [TB] . 



3 Color-ordered amplitudes 

We could now apply this spinor product technology directly to QCD. However, it 
will be useful first to spend a bit of effort analyzing the color structure of QCD am- 
plitudes. It is convenient to divide QCD amplitudes into irreducible, gauge-invariant 
components of definite color structure. We will see that it is most straightforward to 
compute these objects separately and then recombine them to obtain the full QCD 
results. 

3.1 qq-^gg 

Begin with the process qiqu gg- This is similar to the process e~^e~ — ?■ 77 
analyzed in the previous section, except that now there are three diagrams, as shown 
in Fig. |4| including one with a 3-gluon vertex. With the numbering of external states 
as in the figure and all momenta outgoing, the value of the amphtude is 




Hw){-gr''t^) — {ll'^] ■ (e(2) ■ e(3)(2 - 3), 

+eA(3)(2 ■ 3 + 2) ■ e(2) + eA(2)(-2 ■ 2 - 3) • e(3) 



(45) 
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In this formula, t"' and are the color SU{3) representation matrices coupling to 
the gluons 2 and 3, respectively. The third diagram in Fig. |4]has a color structure 
that can be brought into the forms seen in the first two diagrams by writing 

- gf^^H' = {ig) ■ ir^H" = {ig){t^t^ - th") . (46) 

We will find it convenient to rescale the color matrices: T"" = \/2t"', so that the T"" 
are normalized to 

tr [T'^T^] = 6"^ . (47) 



We can then write the amplitude in (45) in the form 

iM = iM(1234) ■ T^T'' + iM(1324) ■ T'T" , (4J 

with 



^M(1234) = (^; 



(l7-e(2)^ ^7 -6(3)4] 

Sl2 



+ — (l7^4][e(2)-e(3)(2-3); 

■S23 



-eA(3)(2 ■ 3 + 2) ■ e(2) + ex{2){-2 • 2 - 3) ■ e(3)] . (49) 



In the second term of (48), M(1324) is given by the same expression with (2,e(2)) 
exchanged with (3,e(3)). 

The elements iM are called color- ordered amplitudes. The complete helicity am- 



plitudes such as (48) are gauge-invariant for any choice of helicities of the external 
particles and for any Yang-Mills gauge group. The two color factors T'^T'^ and T^'T" 
are independent in any non-Abelian gauge group. Thus, the color-ordered amplitudes 
must be separately gauge-invariant. This important observation means that we can 
apply all of the simplifications discussed in the previous section to individual color- 
ordered amplitudes. It will be much simpler to work with these objects rather than 
with the full QCD amplitudes. 

A consequence of this idea is that the individual color-ordered amplitudes should 



obey the Ward identity. It is not difficult to verify this for (49). Replace e(2) by the 
four-vector 2, and use the properties that the external momenta are lightlike, that 
(1 and 4] satisfy the Dirac equation, and that 3 ■ e(3) = 0. Then all terms in the 
resulting expression cancel. 

Other QCD amplitudes can also be reduced to color-ordered structures. Another 
case that will be important for us is the 4-gluon amplitude shown in Fig. |5] This 
amplitude can be written in the form 

iM = iM(1234) ■ tr [T'^T'T'T^] + iM(1243) ■ tr[T"T^T°'T"] 
+iM(1324) ■ tr [T'^T^T^'T'^] + iM(1342) ■ tr [T'^T'^T'^'T^] 
+iM(1423) ■ ti[T''T'^T^T^] + iM(1432) ■ tr[T"T^T"T''] . (50) 
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3 4 



4 1 




+ 

2 1 



Figure 5: Feynman diagrams for gg — )• gg. 




+ f^"f^^^{g^^g^^-g^^g^^) 
+ f^^^f^^^{g^^g^^-g^^g^")] 



Figure 6: The four gluon vertex of QCD. 



I have used the cychc invariance of the trace to move T*^ to the front in each trace. 
Then there are 3! = 6 possible traces, all of which should be included. For a sufficiently 
large SU{N) gauge group, all of these traces are independent; thus, each coefficient 
is a gauge-invariant structure. As in the case of qqgg, these coefficients are given by 
a single function evaluated with different permutations of the external momenta. 

To write the four diagrams in Fig. [5] as a sum of color structures, we need to 
convert color factors in the 3- and 4-gluon vertices into products of matrices. For 



the 3-gluon vertex, this is done through (46), or, by the use of (47), through 



(51) 



For the 4-gluon vertex, we need to apply this decomposition twice. The textbook 
form of the 4-gluon vertex is shown in Fig. |6j Each term can be manipulated as 
follows 
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^^^rparpbrpcrpd rparpbrpdrpC rpbrparpcrpd _j_ rpbrparpdrpc^ (52) 



The full 4-gluon vertex can then be rearranged into 



i|-tr(T"T^T"T'^)[2^^^^^" - g^"" g^'' - g'^'g''^] 



(53) 
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( g^^ (1 -2)^ + g''^ (2-3)^ + g^^ (3-1 )") 



= (2g^^V"-g^V"-g^V') 



Figure 7: Color-ordered Feynman rules for QCD. 



plus 5 more terms corresponding to the other 5 color structures in (50). 



These vertices in the form of traces over T'^'s can be contracted using the identity 
for SU{N) generators 

T^^jTke = 5ii5kj - j^SijSki . (54) 



The coefficients in this equation are determined by the normalization (47), with 5"" = 
A^^ — 1, the number of generators of SU{N), and by the requirement that trT° = 0. 
Using this identity, a product of traces can be transformed into a single trace plus 
smaller factors, for example, 

tr[T''A]tr[T"fi] = tiAB - ^tr[A]tr[5] . (55) 

Using these methods, the value of a Feynman diagram is naturally organized as a 
sum of color-ordered terms. The individual color-ordered amplitudes are computed 
with color- ordered Feynman rules, shown in Fig. [7} 

Using these rules applied to the two diagrams in Fig. |8l it is easy to rederive the 



color-ordered amplitude (49) 



As an example, I will compute the color-ordered amplitude M(1234) that is used 
to build up the four-gluon amplitude. Of the four Feynman diagrams in Fig. |5| only 
the three diagrams shown in Fig. [9] contribute to this color-ordered component. Here 
and in the rest of the lectures, the color-ordered amplitudes in the figures will be 
ordered clockwise. Using the Feynman rules in Fig. [7[ we find 

■'-[e(4) ■ e(l)(4 - 1)^ + e\l){2 1 + 4) ■ e(4) + e^4)(-24 - 1) ■ e(l)] 



^9 s2 



^M(1234)^,^ 
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Figure 8: Feynman diagrams contributing to the first color-ordered amplitude in qq ^ gg. 




Figure 9: Feynman diagrams contributing to the first color-ordered amplitude in gg — )• gg. 



•[e(2) • e(3)(2 - 3)a + eA(3)(2 3 + 2) • e(2) + eA(2)(-2 2 - 3) • e(3)] 



— z 



+— [e(l) ■ 6(2)(1 - 2)^ + 6'^(2)(22 + 1) ■ 6(1) + 6^(l)(-2 1 - 2) ■ e(2)] 
■[e(3) ■ e(4)(3 - 4)a + eA(4)(2 4 + 3) ■ e(3) + eA(3)(-2 3 - 4) ■ e(2)] 
+(-i)[2e(l) • e(3) e(2) • e(4) - 6(1) • 6(2) 6(3) • 6(4) - 6(1) • 6(4) 6(2) • 6(3)]] . (56) 

With some effort, you can show that this expression obeys the Ward identity. That 
demonstrates exphcitly that this color component is independently gauge- invariant, 
as required. 



4 MHV amplitudes 

Now that wc have rcdTiccd the qqgg and gggg QCD amplitudes to managable 
components, it is time to cvalutc these expressions. We will see that the values we 
find fall into surprisingly simple forms. 
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4.1 qqgg amplitude 



Begin witli tlie qqgg amplitude (49). As in the example with e"'"e~77 in Section 2, 
we need to consider all possible cases of gluon helicity. From here on, I will notate 
helicity states by +, — instead of R, L. 

Consider first the case g+g+- Choose the polarization vectors with reference vector 

V2 (12) ' " (13) • ^^^^ 



As with (|4T| ), the first line of (49) is proportional to (11) = 0. Also, in the second 
line of (|49[), each term contains one of the elements e(2) ■ e(3), (I7 ■ e(2)4], and 
(l7-e(3)4], all of which are proportional to (11). So, the entire expression vanishes. A 
similar argument shows that the color-ordered amplitude with two negative helicities 
vanishes. 

The two cases with one positive and one negative gluon helicity are distinct color- 
ordered amplitudes, not related by Bose symmetry. So, both must be computed. I 
will begin with the case g+{2)g^{'i). Choose for the polarization vectors 

' " 72 (12) ' - [43] ■ ^^^^ 

Then all terms in M vanish except for the term in the second line with e(2) ■ e(3). 
This factor is 

.(2) .(3)- 12(13)[42] 

'^^^■'^^^-"27l2)[43f ^ ^ 

Then 

.^^ (1(2-3)4] (13) [42] 
2 (23) [32] ^ (12) [43]^ 

= -.o^^l^fl^ (60) 
(23) [32] (12) ^^^^ 

Again, we can find a simpler form by multiplying top and bottom by (13). Using 
(13) [32] = -(14) [42], we find 

,M=,/, ,f>;f> ^, (61) 

^ (12)(23)(34)(41) ^ ' 

For the case 5'_(2)(7_|_(3), choose for the polarization vectors 
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For this choice of reference vectors, 

e(2) ■ e(3) = , 2 ■ e(3) = 3 ■ e(2) = . (63) 



Then only the first hne of (49) is nonzero. The value of this line is 



^ -W' 1. 2^ (12)[3(l + 2)2)[34] 
2 ^ 2^ (23) [32] (12) [21] 

^ (23) [32] [21] ^ ' 

Now multiply top and bottom by (12) (42). After some rearrangements, all of the 
square brackets cancel, and we find 

'^^-'^ (12) (23) (34) (41) • ^^^^ 

These results are very interesting. The qqgg amplitude vanishes when all of the 
gluon helicities are identical. This followed in a very straightforward way from the 



choice of reference vectors in (57). It is not hard to see that this method extends 
to prove the vanishing of the amplitude for qq plus any number of positive helicity 
gluons. Thus, for any number of gluons, the first nonvanishing tree amplitudes are 
those with one negative helicity gluon and all other gluon helicities positive. These 
amplitudes are called the maximally helicity violating or MHV amplitudes. In the 
above examples, these amplitudes are built only from angle brackets, with no square 
brackets, and have the form 



iM{q_{l)g^{2) ■■■g-{i)--- g+{n - l)g+(n)) = zg 



n-2 (1^) (^^) 



(12)(23) ■■■ ((n - l)n)(nl) ' 

(66) 

where i denotes the gluon with negative helicity. This formula is in fact correct for 
all i, n. The complex conjugates of these amplitudes give the amplitudes for the case 
of one positive helicity gluon and all other gluon helicities negative, 

«M(,.(1),_(2) . . • ■■,-in-m,{n)) = (-^-^ ^^^^ ■ 

(67) 

I will give a proof of these formulae in Section 6. 



4.2 Four-gluon amplitude 

A very similar analysis can be applied to the four-gluon amplitude. Consider first 
the case with all positive helicities. Choose the gluon polarization vectors so that the 
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same reference vector r is used in every case, 
Then, for all i,j, 

e(«)-e(j)~(rr)[jz] = 0. (69) 



By inspection of (56), every term contains at least one factor of e{i) ■ e{j). Thus, the 



entire expression vanishes. 

This argument is easily extended to the case with one negative helicity gluon. The 



amplitude (56) is cyclically symmetric, so we can chose the gluon 1 to have negative 



helicity without loss of generality. Then let 

V2 [21] ' V2 (Ij) ' ^ ^ 

for j = 2,3,4. Again, e{i) ■ e(j) = for all and so the complete amplitude 
vanishes. 

It is not difficult to see that these arguments carry over directly to the n-gluon 
color-ordered amplitudes for any value of n. The tree amplitudes with all positive 
helicities, or with one negative helicity and all of the rest positive, vanish. The 
maximally helicity violating amplitudes are those with two negative and the rest 
positive hehcities. 

It is worth noting that the vanishing of the amplitudes with zero or one negative 
helicity, both for the qq + n gluon case and for the pure gluon amplitudes, is related to 
supersymmetry. QCD is not a supersymmetric theory, but it is an orbifold reduction 
of supersymmetric QCD. In supersymmetric QCD, these scattering amplitudes are 
related by supersymmetry to amplitudes with four external fermions, which must 
contain two negative helicities by helicity conservation along the two fermions lines. 
A precise discussion of these points can be found in [3]. Note that these arguments 
apply only to tree amplitudes. The forbidden amplitudes become nonzero at one loop 
in a nonsupersymmetric theory. 

For the 4-gluon amplitude, all that remains is to compute the color-ordered ampli- 
tude in the case with two negative helicities. By the cyclic invariance of M, there are 
only two cases, that in which the two negative helicities are adjacent and that in which 
they are opposite. As an example of the first case, we can analyze iM(l_2_3+4+). 
Choose the polarization vectors to be 

1 [47'' 1) 1 [47^2) 



V2 [41] ^ ' 72 [42] 

^^3) = ^^ ^^^(4) = ^^ (71) 
^ ^ V2 (13) ^ ^ (14) ^^'^ 
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With this choice, all scalar products of e's are zero except for 



el^j 2 [42] (13) " (13) [42] 



Looking back at (56), we see that the first and third lines are zero. In the second line, 
the only nonzero term is the one that involves e(2) ■ e(3) and no other dot product of 
e's. Thus, 

^M = (^) — (-4)e(2).e(3)2.e(l)3-e(4) 

= (ifM^ (T3) 

^ (34) [41] (41) ^ ' 

Multiplying top and bottom by (12) and rearranging terms in the denominator to 
cancel out the square bracket factors, we find 

iM(l_2_3+4+) = ig \ ,,\,, , . (74) 
^ + +^ y (12) (23) (34) (41) ^ ' 

Similarly, to evaluate iM(l_2_(_3_4+), choose the polarization vectors to be 

1 [47^1) 1 (l7''2] 

V2 [41] ^ ^ V2 (12) 

1 [47^3) 1 (l7''4] 

With this choice, all scalar products of e's are zero except for 

,(9^ (',^ (13) [42] 

^^'^■^^^^ = -(12M ^ ^ 

Again, only the term that involves e(2) ■e(3) and no other dot product of e's is nonzero. 



The value of that term is again given by the first line of (73), which, in this case, 
evaluates to 

- -'^'imSmm ■ '''' 

Multiplying top and bottom by (13)^ and rearranging terms in the denominator to 
cancel out the square brackets, we find 

iM(l_2+3_4+) = tg^^^J^^^^^^ . (78) 
^ + +^ ^ (12) (23) (34) (41) ^ ^ 
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a) 




b) 




Figure 10: Feynman diagrams for quark-quark scattering: (a) ud — )• ud, (b) uu — )• uu. 



The form of ( 74 ) and ( 78 ) strongly suggests that the general form of an ra-gluon 

(79) 



MHV amplitude is 



(12)(23)---((n-l)n)(nl) 
The corresponding formula, exchanging positive and negative helicities, is 



iM{g_{l) ■■■g+ii)--- g^j) ■ ■ ■ g-{n)) = (-1)^^?""' 



[12][23]---[(n-l)n][nl] ' 



^0) 



I will give a proof of these formulae for all j, n in Section 6. The formula (79) was 



discovered in 1986 by Parke and Taylor [17j . This was the original breakthrough that 
gave the impetus for all of the work discussed in the latter sections of this review. 



5 Parton-parton scattering 

Before going deeper into the theory of the MHV formulae presented in the previ- 
ous section, I will present a simple application of these formulae. The basic ingredient 
for collider physics is the set of tree-level cross sections for parton-parton scattering. 
These are straightforward to derive from the QCD Feynman rules, and yet the cal- 
culations can be tedious for students. In a one-year course in quantum field theory, 
this subject generally arises just at that point in the year when the professor's family 
needs a weeklong ski vacation. Thus, in the textbooks, the derivation of these formu- 
lae is typically left as an exercise for the students without detailed explanation. In 
this section, I will show that the MHV formulae make these derivations trivial. 



5.1 Four-fermion processes 

Begin with quark-quark scattering. For scattering of quarks of different flavor, 
there is only one Feynman diagram, shown in Fig. ITota). The value of this diagram 
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can be found from (26). In this formula, all momenta are outgoing, so we simply need 
to cross the correct particles into the initial state. Which particles should be crossed 
depends on which helicity amplitude we are considering, for example, uidL u^di 
or uidfi — )■ uidfi. We also need to include the QCD color matrices. With the 
numbering of lines as in the figure, with the arrows indicating the direction of left- 
handed fermions, the matrix elements are 

iM{uLdL ^ ULd^) = 'a'^^l^f^Tl, 

(81) 

The scattering of identical quarks, for example, ulUl — )■ ulUl, is given by two 
Feynman diagrams, shown in Fig. [Io|(b). For this case 



iM{uLUL ulUl) = ig 



(34)^ (43)^ 

\ / rpa rpa \ / rpa rpa 

1MC2^ ' ' (14)(23) ^1 



^2) 



The extra minus sign in the second term comes from interchange of identical fermions. 
To compute cross sections, we square these expressions, sum over final helicities and 
colors, and average over initial helicities and colors. The interference term in the 
square of (82) requires 

(13)(24)((14)(23))* = (13) [32] (24) [41] = -(13) [31] (14) [41] = -tu . (83) 

By parity, the cross sections are unchanged when all helicities are reversed. The 
cross sections for antiquarks can be computed by applying crossing symmetry, so the 
formulae above are all that we need to cover all of the possible cases. 

For the color sums and averages, here and in the later calculations in this section, 
we will need the formula T'^T'^ ~ traces 

tr [T'^T'^T^T''] = — , tr [T^T^'T'^T*] = - ^ . (84) 
3 3 



These are easily proved using ( 54 ) . The second trace appears in the interference term 
in the square of (82). For a general SU{N) gauge group, these results are 

tr [T^T'^T'T^] = ^ — , tr [T'^T^T'^T^] = ^ ^ , (85) 

so that the second color sum is suppressed by l/N"^. This is an example of the 
general result that interference terms between different color structure are suppressed 
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Figure 11: Notation for the process ulur QigR- 

by l/N"^. That result follows from the broad picture of counting presented by 't 
Hooft in [18]. 

Assembling the pieces, we find 



da 



dcos9 
da 

dcos9 



{ud —7- ud) 



[uu — )■ uu) 



2-nat 



9 s 



+ 



9 s 



— + — 



2s2 
3 tu 



(86) 



where 6 is the scattering angle in the parton-parton center of mass frame. The second 
of these formulae is to be integrated over cos6 > only. Crossing these formulae into 
other channels, 



da -J -J 2 7ra^ 
{ud ud) = - — - 
9 s 



dcos 9 

da 
dcos 9 



[uu — )■ uu) 



+ u^' 
2-nal \s'^ + u^ 



9 s 



+ 



V + u^ 



2_uh 
3 St 



(87) 



5.2 qqgg processes 

Next, consider qq — t- gg. The q and q necessarily have opposite helicity. The 
gluons also have opposite helicity. The two cases of gluon helicity are related by 
gluon interchange or t ^ u. So it suffices to consider the case shown in Fig. [TT} 
ul{4:)ur{1) — >■ gL{2)gfi(3). For this case, the color-ordered amplitudes are MHV 
amplitudes, and so our results from Section 4 give 

(13)^(43) , (13)=^(43) , 1 

(12) (23) (34) (41) (13) (32) (24) (41) J ' ^ ' 



iM = ig"^ 
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The square of this amphtude, traced over color, is 

tr[M|' 



9 
64 



64 ,tu f 
9t2 



2^ 



3 M 4 s'^ 



^9) 



From the first hne to the second, I have used w^+t^ = {u+ty —2ut = s'^—2ut. Adding 
the case with the two gluons interchanged, and including the factors for averaging 
over initial colors and helicities, we find as the final result, 



da 



dcosO 



[uu gg) 



16 Tia 
27" 



2 r 



t U 
U t 



4 s2 



(90) 



This formula should be integrated over cos^^ > only. Crossing this formula into 
other channels gives 



da 
dcos9 

da 
dcos6 



{ug ug) 
{.99 -> uu) 



9^ 
1 vra 
12' 



s u 9 s + n 

+ 7 ,o 

u s A 



2 r 



t U 
U t 



9t2+^ 
4 s2 



(91) 



5.3 Four-gluon processes 

Finally, we need to compute the gg — )■ gg scattering amplitude. This can be done 
directly from the formalism we have already developed, but it is useful to add one 
additional trick. 

Extend the gauge group from SU{N) to U{N) by adding an extra generator 
T° = \/\/N. Since the 3- and 4-gluon vertices are proportional to the structure 
constants /"^^ and the extra U{1) generator commutes with the SU{N) generators, 
the additional boson gives no contribution to the cross sections. However, adding this 



term simplifies the color algebra. The contraction identity (54) now takes the simpler 
form 

T^^jTke = SiiSkj . (92) 

We can use this identity to compute the square of the 4-gluon amplitude ( [50) ). The 
squares of the individual color traces are equal to 

tr[T'^T*T"T^] tr [T'^T^T'T"] = A^^ = 81 . (93) 

The cross terms are proportional to 
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= 



Figure 12: The U{1) Ward identity for color-ordered amplitudes. 



Then, letting I, J = 1, . . . , 6 index the six possible cyclic orderings of the gluons, we 
find 



tr[|7W|2] = 81 5] |M(/)p + 9 ^ M{I)M{jy 
= 72 5]|M(J)p + 9|^M(/)p . 



(95) 



I explained above that, when we add the extra generator to extend SU{N) to 
U{N), the value of tr|A^p cannot change. The vanishing of the coupling to the r° 
gauge bosons is expressed in the color-ordered amplitudes as the Ward identity shown 



in Fig. 12 The sum over all possible orderings of a U{1) boson in a color-ordered 
n-gluon amplitudes must vanish. An example of this identity, for 4-gluon amplitudes, 
is 

iM(0_l_2+3+) + iM(l_0_2+3+) + zM(l_2+0_3+) = (96) 
We can easily check this from the explicit expressions for MHV amplitudes. Evalu- 



ating the left-hand side of (|96|) gives 



1 



(01)(12)(23)(30) (10)(02)(23)(31) (12)(20)(03)(31) 
z^?2(oi)^[<02)(13) + (03)(21) + (01)(32)] 



(97) 



(01)(12)(23)(30)(02)(13) 
The expression in brackets in the second line vanishes by the Schouten identity (19) 
The U{1) Ward identity implies that 



Then (95) becomes 



EM(/) = o. 



tr[|A^r] = 72E|M(/)P 
I 



{91 



(99) 
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Now we must evaluate this expression for all possible choices of the external parti- 
cle helicities. Of the 16 possibilities, only 6 involve 2 positive and 2 negative helicities. 
All other choices give zero. For the color-ordering 1234, the nonzero amplitudes are 

,4 



|Mfl_2_3+4, 



|Mfl_2,3_4, 



9 



(12) (23) (34) (41) 
(13)^ 



9 



(12) (23) (34) (41) 
(14)^ 



I 



|M(l_2+3+4_)P = , , , \ ' , , = , (100) 
' ^ + + ^1 y (12) (23) (34) (41) ^ s2 ' ^ ^ 

and the parity conjugates of these amplitudes. The sum over color orderings sums 
over these results crossed into all possible channels. Thus 



EEtr|M(I) 

h I 



16/ [3 - 



t2 

su 



+ 



+ 



u 



+ 



+ 



ut 

T2 



st 



Including the factors for initial-state color and helicity averaging, we find 



da 



-X99 99) 



9 vra^ su ut st , 
4~^ " ~ ^ ~ 



(101) 



(102) 



dcosO 

This formula should be integrated over cos 9 > only. This completes the calculation 
of all of the tree-level 2 — )■ 2 parton-parton scattering cross sections. 



6 Britto-Cachazo-Feng-Witten recursion 

Now I will develop some methods that will allow us to prove the MHV formulae, 
and also to use these results to compute the more compex non-MHV amplitudes. The 
general approach will follow the idea of Cachazo, Svrcek, and Witten P^PU] that we 
should consider the tree-level color-ordered amplitudes as analytic functions of the 
angle brackets and square brackets. We can then analytically continue scattering 
amplitudes to complex momenta while keeping the property that all external legs 
are on mass shell. At the end of the calculation, to extract physical results, we will 
specialize to the values such that {ij) = {[ji])*. 

6.1 Three-point scattering amplitudes 

This freedom to discuss QCD scattering amplitudes for arbitrary complex mo- 
menta allows us to write on-shell three-point scattering amplitudes. For three-point 
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1 



Figure 13: Diagrams contributing to the three-point qgq and ggg amphtudes. 



amplitudes, shown in Fig. 13, momentum conservation implies that 1 + 2 + 3 = 0. 
Then 

= l2 = (2 + 3)2 = 22-3 , (103) 

and similarly the other two Lorentz products vanish, so it would seem that there are 
no invariants for these amplitudes to depend on. However, in terms of square brackets 
and angle brackets, 

22 - 3 = S23 = (23) [32] . (104) 



Thus, with complex momenta, we can satisfy (103) by having [32] = while keeping 
(23) nonzero. 

With this idea in mind, I will compute the color-ordered amplitude for gluon 
emission from qq, i'M.{q_{l)g_{2)q_^_{3)). The amplitude comes from one Feynman 
diagram, shown in Fig. [T3j The value of this diagram is 

with r an arbitrary lightlike reference vector. Applying the Fierz identity, and then 
multiplying top and bottom by (12), this rearranges to 

(12)[r3](21) (12)' , , 

^ [r2](21) ^ (31) ^ ' 

The reference vector r cancels out, and the final result is just the MHV amplitude 
with n = 3, 

M(,_(l)9_(2)5,(3)) = .9^|g|y. (107) 

A similar calculation gives the form of the three-gluon vertex. 
tM{g_{l)g^i2)g+{3)) 
= ^[e(l)-e(2)(l-2)-e(3) 

+e(2) • 6(3)(2 - 3) ■ e(l) + 6(3) ■ e(l)(3 - 1) ■ 6(2)] . (108) 



28 



Choose for the polarization vectors 

eni)---^ e-m--^!^ e^(3)-^^-^ (109) 



so that e(l) ■ e(2) = 0. Then (108) evaluates to 

ig 1 (0+[r3](s2)[r(2-3)l) + [r3](sl)[r(3-l)2)) 
V272 [rl][r2](s3) 

[rl\ [r2\ 

If we multiply top and bottom by (12)^ and rearrange the denominator, the factors 
of (r3) cancel, and we find 

(12^^ 

''^to-<^)«-P'«^P)'-'« (12)W)(31) - 

For reference, the corresponding formulae for two positive and one negative helicity 

are 

M(,_(1K(2)?,(3)) = -Jpllgj , tMMl)u(2)9-{3)) = ■ 

(112) 

I note again that all other n-gluon and qq+n gluon amplitudes with only one negative 
or positive helicity are zero. Only these cases are nonzero, with the MHV values. 



6.2 BCFW recursion formula 



It would be wonderful if we could use these three-point functions as building blocks 
for the construction of general tree amplitudes. The usual understanding is that we 
need off-shell three-point functions to build up general amplitudes. However, this 
common-sense idea is evaded by a beautiful strategy of Britto, Cachazo, Feng, and 
Witten EH. 



Consider a color-ordered amplitude iM(l ■ ■ ■ n), as shown in Fig. 14 Choose two 
legs i, j, and choose a value of z, a complex variable. Now define new momenta i and 
j by shifting the square bracket of i and the angle bracket of j, 

i) =i) i]=i] + zj] 

J)=j)-z^) J]=J], (113) 
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Figure 14: Notation for the discussion of the Britto-Cachazo-Feng recursion formula. 

The momenta i and j are given by 

i = i)[i + zi) [j j = j) [j - zi) [j . (114) 

Notice that i + j = i + j, so momentum conservation is respected. For general z, i 
and j have unphysical complex values. However, since is proportional to {ii) = 0, 
and similarly, = 0, the new momenta remain on shell. 

As a rule of thumb, one should shift the square bracket of an external state with 
negative helicity and the angle bracket of an external state with positive helicity, that 
is, = (— , +). This minimizes the number of factors of z in the numerator. I will 
give a more careful analysis of this point below. 

Call the amplitude evaluated at the modified momenta iM(2;). Now consider the 
integral 

^ tM{z) (115) 



27ri z 

taken around a large circle in the complex z plane. If ilSA{z) — > as 2; — > 00, then this 
integral vanishes. In that case, we can set to zero the value obtained by contracting 
the contour and summing over the poles that it encloses. There is an obvious pole at 
z = 0. The residue of this pole is exactly 

iM(^ = 0), (116) 

the amplitude that we wish to evaluate. 

The other poles come from the denominator of the color-ordered amplitude. We 
are discussing tree amplitudes, and so these poles will come from the propagator 



involved in a factorization of the amplitude, as shown in Fig. [15] The denominator 
of this propagator is 



b 

\2 



= Sa...biz) = kmY , (117) 
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Figure 15: A factorization of the amplitude in Fig. 14 that gives a pole in the amplitude 
iM{z). 



the invariant sum of momenta on the right, including any terms from the shifts (114). 
If both i and j are on the same side in the factorization, Sa--b is idependent of z and 
there is no pole. If i is on the left side in the factorization and j is on the right, 

b 

Sa-biz) = ( J2 ^rn- zi)[j f 
m=a 

b 

= Sa...b{0) - z{t{J2 k^)j] + z'mjj] . (118) 

a 

The last term is zero. The denominator then has a simple pole at 



(119) 



The residue of this pole is 
1 



iM{{b + 1) (a - l)Q) ,., iM{Qa ... j ... 6) 



-iM((6 + 1) (a - l)Q)—iM{Qa ... j ... 6) (120) 

Sa-b 



The amplitudes in each term are evaluated with the value z = z*, given by (119), 
appropriate to the cut that is chosen. 

The result of this analysis is that, if we have chosen the shifted momenta i and j 
so that iM.{z) vanishes at infinity, then 

iM{l...n) =^iM((6 + l)...?...(a-l)(-g)) — iM{Qa . . .j . . .b) . (121) 

a,b ^"■■■■^ 

where Q = -El ^m- This is the BCFW recursion formula. An arbitrary color- 
ordered helicity amplitude can then be evaluated by breaking it down into amplitudes 
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with a smaller number of legs. The amplitudes on the right-hand side of (121) are 
evaluated with on-shell but complex-valued momenta. 

The recursion must be run until we can immediately evaluate the amplitudes on 
the right-hand side. Since all nonzero amplitudes with five external legs or fewer are 
either MHV or the conjugates of MHV amplitudes, only a few steps of the recursion 
are needed in practice. 

For QCD with all massless particles, the complete solution of the recursion is 
known. However, the result is not simple, so, in keeping with the philosophy of these 
lectures, I will refer you to the reference P^23f24j 

There is one subtlety that should be clarified to evaluate the right-hand side of the 
BCFW recursion formula. The amplitudes involve angle brackets and square brackets 
of the complex momentum Q. In the examples given later in our discussion, I will 
evaluate these brackets by assembling them into complete factors of the momentum 
Q. To do this, we will need to relate the brackets {—Q)) and {—Q)] in the amplitude 
on the left to Q) and Q\. It is consistent always to take 

(_Q)) = ^Q) (-g)] = ^Q] . (122) 

One special circumstance should be noted. If the line on which the amplitude factor- 
izes is a fermion propagator, the value of this propagator is 

Then one of the brackets in the left-hand amplitude is a Q], not a {—Q)]- To com- 
pensate for this, we need to add a factor {—i) for a cut through a fermion propagator. 

Finally, we need to discuss under what circumstances the amplitude zM(z) will 



vanish at infinity. I claim that, for the shift given by (113) in which i is shifted in 
the square bracket and j is shifted in the angle bracket, the helicity choices = 
(—, -|-), (—,—), -|-) give the BCFW recursion formula, while in the case = 
— ) there are extra terms from the integral at infinity. A transparent argument 
for this claim has been given by Arkani-Hamed and Kaplan |25j . 

In the limit 2; — )■ 00, the dominant momentum flowing through the diagram is 
that induced by the shift. Define 

q = t)[j or q^ = {tYj] (124) 

Then the large momentum is zq. Note that = 0. It is useful to picture the diagrams 
as containing one line whose propagators carry momenta zq + k close to this large 
momenta. This line attaches to propagators that carry momenta of order 1. This 



32 



structure is illustrated in Fig. 16 Each propagator has a denominator of order 

{zq + kf = 2zq-k + k^ = 0{z) . (125) 



Each vertex is at most linear in momentum, and so is at most of order z. Then 
each combination of vertex and propagator is of order z^. The chain of vertices and 
propagators has one extra vertex, so this is of order z. 

Now consider the polarization vectors. We can choose i or j as the reference 
vector, whichever gives a nonsingular result. Then 

. _ z{jYj] + Ul^i] , _ ^ [^^^^) 

For the choice (i, j) = (— , +), from the 2;-dependence of the polarization vectors and 
the internal vertices and propagators, we see that the amplitude cannot be larger 
than 0{1/ z). For the choice (i, j) = (—,—), there is apparently a term of order z^. 
However, it is not possible to build this term, since 

e_«-e_(j) = g.e_« = g-e,(j) = 0. (127) 

The first nonzero term is one in which both e_(i) and e_(j) are dotted with 0{1) 
vectors instead of zq^, and this term is down by two powers of z, giving in all an 
amplitude of of 0{l/z). An analogous argument holds for the (+, +) case. The 
amplitude in the (+, — ) case is irremediably 0{z'^). 

A parallel argument can be carried out for the case in which the shifted lines are 
fermions. In the fermion case, the vertex involves no momenta and is of order z", but 
the propagator is also 0{z^). The leading term in the propagator is proportional to 
^, which annihilates both i) and j]. Then, again, only the (+, — ) case can be nonzero 
as 2; — )■ 00. There is one exception to this rule: If we shift on momenta i and j at 
opposite ends of the same fermion line, the diagrams with only a single vertex and 
no internal propagators on this line can give a term of 0{z^). 

The paper taking a somewhat more sophisticated approach, generalizes this 
conclusion to other theories, including gravity. 



6.3 Proof of the MHV formula 



As a first application of the BCFW recursion formula, I will give a proof of the 



MHV formula (79) for n-gluon amplitudes, using an argument originally given by 
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Figure 16: Analysis of the z ^ oo behavior of a multigluon tree diagram. 




+ 




+ ... + 



(N-l)+ (k+l)+ k+ ■■■ 'j^ 




Figure 17: BCFW cuts contributing to the evaluation of the n-gluon MHV amplitude. 



Risager 
method. 



The other MHV formulae in Section 4 can be proved by the same 



The proof will proceed by induction. We have already verified this formula for the 
n = 3, 4 MHV gluon amplitudes. So, I will assume that the MHV formula is correct 
for the case n = N — 1 and use that hypothesis to evaluate the n = N gluon MHV 
amplitude. 

The iV-gluon amplitudes are cyclically invariant. So, without loss of generality, 
choose the gluon 1 to be one of the two gluons with negative helicity. Define the 
z-dependent amplitude in the BCFW procedure by the shift 



i] = 1] + ^2] 2) = 2) . 



(128) 



The BCFW recursion formula then evaluates the A^-point MHV amplitude in terms 
of factorized diagrams with 1 on one side of the factorization and 2 on the other. 



as shown in Fig. [T7| In each diagram, we must assign a helicity to the intermediate 
gluon. There are two choices: + outgoing from one side and — outgoing from the 
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Figure 18: The three-point amphtude from the first diagram in Fig. 17 



other side, or vice versa. In most cases, both choices give zero, since one side or the 
other will be an n-gluon amplitude with n > 3 and one or zero — helicity, which 
vanishes by the arguments given in Section 4. The two diagrams that are potentially 
nonzero are those in which one side or the other is an amplitude with n = 3. These 



are the first and last diagrams in the series shown in Fig. 17 



The first diagram on the right-hand side in the figure involves the 3-point vertex 



shown in Fig. 18 The denominator of the exposed propagator is 

= {lN)[N{l + z2)] (129) 
However, the 3-point vertex has the form 

- 4 

-tg^rJ^ ^ = +tg\N(l + z2)] . (130) 

Then vertex cancels the pole in z from the propagator, so the diagram actually has 
no pole. 

The only nonzero contribution to the sum then comes from the last diagram. Its 
value is 

,,--3 (liV(-l) - ( [23]^ ^31) 

(lQ)(g4)(45)---(7Vl)(23)[32]^ [g2][23][3g] 

with 

g = -2)[2-3)[3 + 2l)[2 . (132) 



From (132), we derive 

(ig)[g3] = -(i2)[23] 

(4g)[g2] = -(43)[32] . (133) 
The factors of [23] all cancel, and we end up with the result 

ia^-^ fl34l 

^ (12)(23)(34)(45)---(Ari) ^ ' 
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which is exactly the Parke- Taylor amplitude for the case ofn = N legs. By induction, 
this formula applies for all n. 

The MHV formula for n-gluon amplitudes with 2 + helicities, and the MHV 
formulae for amplitudes with 2 fermions and {n — 2) gluons can be proved by following 
the same strategy. 



7 pp ^ + partons 

As an illustration of this technology, I will now derive the complete set of formulae 
needed to compute the cross section for vector boson production with up to 2 partons 
at a hadron collider. The methods I have discussed make the computation of the 
cross section for a vector boson plus one parton truly trivial. Only some particular 2 
parton amplitudes will require a little work. 

For definiteness, I will write formulae for production from u and d quarks. 
The contributions from other light flavors, and the corresponding formulae for W~ 
and production, can be derived from these by small modifications of the prefactors. 



7.1 ud^W+ 

The coupling of the to quarks and leptons is 

SL = ^W7{dYPLU + h^PL^) + h.c. , (135) 
V2 

where Pl = {1 — 7^)/2 is the left-handed projector, and Qw is the weak interaction 
coupling, satisfying 

The W~ field appears because this is the field that creates the . This interaction 



leads to an amplitude given by the Feynman diagram shown in Fig. 19 



iMiud vt) = i —u(l)a''u(2)u(3)w.u(A) 

= '9l- . ^ (31) [24] . (137) 

We can treat this formula in one of two ways. One way is to integrate over the full 
phase space of the final-state leptons. This will give a broad mass distribution for the 
leptons, on top of which the will appear as a resonance. 
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Alternatively, since the is a narrow resonance, it is typically a good approx- 
imation to compute amplitudes for a real boson on its mass shell. In this case, 
you might think it is easiest to sum over the polarization vectors. However, we 
will obtain simpler formulae if we retain the final-state lepton spinors. As a bonus, 
retaining the spinors will preserve the information on the polarization. This is 
valuable, first, because the polarization gives useful information on the dynamics 
of the W'^ production and, second, because the experimental acceptance for a 
depends strongly on its polarization. 

The matrix element for the coupling to leptons is proportional to the current 
u{i')a^v{i~^). This matrix element must be squared and integrated over the direction 
of the leptons in the rest frame. That integral is 

1"' = l^{ir2]{2Yl] ■ (138) 

To evaluate this integral, note that, if q'^ is the momentum, q = 1 + 2, then 
q^Ii^'^ = qji^'^ = 0. Further, 

= / ^ 2(12) [12] = -2g2 = -2ml. ■ (139) 

From these requirements. 



I^'^ = - ml^ig^^ - ^-^) . (140) 



The quantity in parentheses in this equation is the usual sum over on-shell 
polarization vectors. Thus we can represent this sum as 



i:<MCW = ^/S(10-2](2yi], (141) 



This gives a simple procedure for computing cross sections with a final on-shell 
p3]: (1) Write the Feynman diagrams with an internal propagator and the 
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final state z/(l)£+(2). (2) Remove the factor (gw/\/2) /(si2 — myy + iniw^w) and 
put tlie momentum (1+2) on sliell. (3) Evaluate the amplitude with spinor 
products. (4) Square and integrate over phase space including the on-shell , and 
add an integral over the lepton direction in the rest frame: 



3 r dVt 



dn„,„../dn„«^y-. (i42) 



The cross section for W'^ production is typically already a multi-dimensional Monte 
Carlo integral, so the computational price of adding two additional, well-behaved 
integrals is small. 

To illustrate this formalism, I will work out the cross section for ud — > . I use 
the matrix element (137) in the form 

iMiW^du) = ^(l7'^2](37'^4] = V2g^{31)[2A] . (143) 

The phase space integral / c/IIo+vk contains 1-body phase space 

J dUi = 27rS{s - m^) . (144) 
Then, averaging over initial colors and spins, we find 



1 Sna,,, _ o , 3 f dVt 



^ . 2rf(. - ™'^) ■ ^ y ^ . (145) 



The integral over dVt implements the familiar (1 — cos^^)^ angular distribution of the 
decay lepton with respect to the u quark direction, signalling a with left-handed 
polarization. Evaluating the integral, we find 

a{ud^W+) = ^^5{s-ml) , (146) 
the familiar expression for the Drell-Yan cross section. 



7.2 MHV amplitudes 

To use this formalism for computing production cross sections, we need to 
be able to compute the amplitudes. We might hope that many of these amplitudes 
belong to an MHV series and are therefore trivial to evaluate. In fact, it is so. 

To explain this, I will first compute some amplitudes for the related process 



e^e qq + ng. We have seen the first of these, for zero gluons, already in (26) and 
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Figure 20: Notation for the amplitudes for eje^ —^qL + ng + qj^. 



(28). This amphtude seems to belong to both an MHV and an anti-MHV series. For 
convenience, I will renumber the external legs to correspond better to the formalism 
that we used in Section 4. On the electron side, I will take the outgoing left-handed 
line (the incoming e^) to be line 1; on the quark side, I will take the outgoing left- 
handed line (the outgoing q^) to be line 3. I will notate a general amplitude for 
^R^L ~^ iL + ng + qj^ in the same way, as shown in Fig. 20 Then the n = case gives 

^M = -2^e^^^ = -2^e^^ (147) 
'■^^ ^ (12)(34) [12][34] ■ 

The electric charges gives an additional factor {—Qf); I omit this here and in the 
following. 

It is not difficult to work out the next case, n = 5 or 1 gluon, using the methods 
described in this review. In the case of a + helicity gluon, it is easiest to take the 
reference vector for the gluon (line 4) to be 3; for a — helicity gluon, take the reference 
vector to be 4. Then one finds 

/13\2 

tMie4l)e42)q43)g4A)q45)) = -2ie^gs ^ ' 



'(12) (34) (45) 
[241 2 

zA^(e_(l)e+(2)g_(3)(7_(4)g+(5)) = 2ze2^/.pi^^ , (148) 

times the factor (— Q/) and the color matrix T'*. The sum of the squares of these 
expressions gives 

Y.\M\' = AeVs'^^^ . (149) 
This reaction is usually described using the kinematic variables 

_ 2^3 ■ q _ 2^4 ■ q _ _ 2k^ ■ q 

q^ q^ q^ 
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where g = 3 + 4 + 5 = —(1 + 2). Each Xi is the ratio of that particle's momentum in 
the center of mass frame to its maximum value Also, 

■S34 = g^(l - a;?) , S45 = q'^{l- Xq) , S35 = g^(l - Xg) . (151) 

The standard expression for the cross section for gluon production in e~^e~ is 



da _ 47r«2 2a, x^ + x| 



doc Q doc q 



3s 



37r 



3'(jr)(l Xq^ 



(152) 



The last factor is almost exactly of the form of ( 149[ ) after substituting (150), (151). 
By reintroducing Q/, summing over final colors, averaging over helicities, and aver- 
aging over the relative orientation of the final state plane with respect to the beam 
axis, it is not difficult to complete the derivation of (152) from (149). 



It is tempting to guess that the formulae (148) generalize to the color-ordered 
amplitudes for all + helicity and all — helicity gluons according to 

2A^(e_(l)e+(2)g_(3)(?+(4) ■ ■ ■ g+{n - l)g+(n)) 



(12)(34)(45)---((n- l)n) 
zA^(e_(l)e+(2)g_(3)(7_(4) ■■■g-{n- l)q+{n)) 

[2nf 



n+l 



[12] [34] [45] ■ ■ ■ [(n - l)n] ' 



(153) 



This is in fact correct. These new MHV formulae can be proved using the method 
of Section 6.3. 

Finally, we need to make the connection between the QED amplitudes given here 
and the amplitudes that we need for production. The only change needed is to 
substitute the photon propagator with a propagator by multiplying by the factor 



9n 



S12 



2e2s 



12 



rriL + imw'^w 



(154) 



and then to take the on-shell according to the prescription given below (141). I 
will write the on-shell amplitudes in the form 



i7W(Vr+(12);3' 



■n] 



V2gu,9s-^T''' . . . T"" ■ iM(12; 34 ■ ■ ■ (n + 4)) 

+other color structures (155) 



Then for the cases discussed above, 



M(12:3 - ■ -n) 



[12](13)V(12)(34)(45)---((n-l)n) 
(-l)"(12)[2n]V[12][34][45] ■■■[{n- l)n] 



MHV 
anti-MHV 



(156) 
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7.3 W+ + 1 parton 

The amplitude for ud — + g+ contains only one color factor. The correspond- 
ing color-ordered amplitude, following from (148) or (156), is 

(34) (45) • 

Using this result, it is straightforward to compute the lowest-order cross section for 
production of + 1 parton. There are three processes that must be considered: 

ud W^g , ug W^d , gd W^u . (158) 



For all three processes, the amplitude is an appropriate crossing of that given in (157). 

For ud — )■ W^g, after averaging over initial colors and helicities and summing over 
final colors and helicities, we find 

a{ud W+g) = — J dUi+w 

|M(12;rf(3)(7+(4)n(5)|2 + |M(12;d(3)(?_(4)^(5)|2l . (159) 



I have used the color sum 

tilT^T^] = 8 . (160) 
If 0* is the polar angle in the W~^g center of mass system. 



" I d cos 9^ f dn 



with s = S5, and 



2 '2 

|M(12;d(3)(7+(4)lZ+(5)r + |M(12;rf_(3)(7_(4)^Z+(5)r = '-^^^^ . (162) 

-534-345 



Then 



a{ud ^ W^g) = / / ^(1 - "^) r'M±4, ) . (163) 



3s J 2 J 47r s V SdgSgu 

Similarly, 

a{ug WU) = a{g-d W^u) = ^ / fs^+^A 

As J 2 J An s \ Sdg 



vra^as f dcos9^ f dfl ^w^^s'^d + ^l^ 

(164) 

where the Sij invariants are evaluated with appropriately crossed momenta. In the 
two cases, respectively, Sgu = s and Sdg = s. 
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7.4 Cross section formulae for + 2 partons 



For the case of W + 2 partons, many reactions must be included, and each of 
these has a nontrivial color structure. The reactions are all of the form (2 partons) 
— >■ (2 partons) and so they are in 1-to-l correspondence with the 2-to-2 parton 

processes that we considered in Section 5. Further, since the is a color singlet, 
the color structure and color sums are also just those that we met in Section 5, and 
we can borrow the color algebra from that discussion. This will allow us to write the 



cross sections in terms of the matrix elements M normalized as in (155) 



We can begin with 4-fermion reactions. For the simplest case of a non-identical 
spectator quark, the cross section is given in terms of the amplitude M(12; 34; 56) 
containing three fermion lines: the leptons (12), the left-handed quark emitting the 
(34), and a second left-handed quark (56). This cross section is 

a{us W+ds) = ■ 1287r=^ J dU2+w ||M(12; 34; 56)^ + |M(12; 34; 65)|H , 

(165) 

The sum of matrix elements takes care of the polarization sum for the spectator 
quark. For spectators identical to the final quark. 



aiud ^ W+dd) = ■ 1287r=^ / dU^+w 

9 s J 



|M(12; 34; 56)1^ - jRe[M(12; 34; 56)*M(12; 54; 36)] 



|M(12;34;65)|^]| . (166) 



I have adjusted the normalization so that the integral can be taken over all of phase 
space. The cross section a{uu — )■ W^du) involves a similar sum over matrix elements. 
Cross sections involving antiquarks have the same form, except that the amplitudes 
M(12; 34; 56) should be crossed appropriately. 



The cross section for the reaction ud — )■ + 2g is given, analogously to (90), by 

„2 



a{ud ^ W^gg) = ■ 128vr3 / dll^^w E 



|M(12; 34i5j6)p - ^Re[M(12; 34i5j6)*M(12; 35^4^6)]) , (167) 
8 J 

where i,j = +, — are the helicity states of the final gluons. The cross section is 
again normalized to be integrated over all of phase space. The cross sections with 
initial-state gluons are given by 

1 Q^-) 11 
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||M(12; 34i5j6)|2 + |M(12; 35^4^6)1^ - ^Re[M(12; 34i5j6)*M(12; 35^-4^6)] | 
a{gg ^ W^du) = ■ ISStt^ / rfH^+^y ^ 

||M(12; 34i5j6)|2 + |M(12; 35^4^6)1^ - ^Re[M(12; 34i5j6)*M(12; 35^-4^6)] | , (168) 

where the matrix elements M(12; 3456) are crossed appropriately. 

There is no analogue of the gg — ?■ gg process, since at tree level the does not 
couple directly to gluon lines. So now we have all of the cross section formulae, and 
it only remains to compute the amplitudes. 



7.5 Amplitudes for + 2 partons 

The cross section formulae in the previous section involve two sets of amplitudes, 
those for Wqggq and those for Wqqqq. These must be evaluated for all cases of the 
final-state helicity. However, each quark line must have a left-handed fermion at one 
end and a right-handed fermion on the other end, so the cases that must be considered 
are only those of different gluon helicities. Of these, two are MHV or anti-MHV, so 
only two amplitudes with gluons, plus the four-fermion case, need to be computed. 
All of those computations are easily done by the BCFW method. 

The two amplitudes that are trivially known are: 

M(12;ci_ (3)^44)^45)^46)) = ^^^^ 

M(12; d.{3)g.{A)g_{5)uAQ)) = • (169) 



I will now take up the case of the amplitude i'M{12] d^{3)g+{4:)g^{5)u+{6)) . To 
evaluate this amplitude, carry out a BCFW shift on the legs 5 (with square brackets) 
and 4 (with angle brackets, according to 

5] = 5] + ^ 4] , 4) = 4) - z 5) . (170) 



Then the BCFW recursion formula gives the identity shown in Fig. 21 The vertical 
dashed lines show the two BCFW cuts. 

I will evaluate the first of these cuts explicitly. The amplitude to the left is an 
MHV amplitude for Wugd. The amplitude to the right is a nonzero MHV 3-point 
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+ 




Figure 21: Evaluation of the amplitude zM(12; d-{3)g^{4:)g^(5)u+{6)) by BCFW recursion. 



vertex. Together, these give 

[12](13)^ 



^ . (Q5) (65) 



(34)(4(-g)) (g5)(56)(6g) 



(171) 



The factor (— i) is associated with the fermion cut; see below (123). Using 5) = 5) 
and —Q) = iQ), this simphfies to 



[12](13)^(5g)' 

(34)s56(4g)(6g) 



(172) 



The momentum Q takes the value 

g = g)[g = -5)[5-6)[6-^5)[4. (173) 

This should be evaluated at the value z = z^, for which = 0. That is, 

S56 [65] 



(5(5 + 6)4] [64] ■ 



(174) 



" 2 " 2 

To evaluate the angle brackets of g, multiply (172) by [g4] /[g4] and use (173) to 



evaluate Q)[Q. This gives the identities 



Also, using (174) 



(5g)[g4] = -(56)[64] 
(6g)[g4] = -(65)[54] 

(4g)[g4] = -(S45 + S46 + Sse) = -S456 

(34) = (3(4 + 5)6]/[46] . 



(175) 
(176) 
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With all of these simplifications, the expression for the cut becomes 

[12](13)'[46]' 
[45][56]s456(3(4 + 5)6] " 



(177) 



After evaluating the second cut in the same way, we arrive at the following ex- 
pression for the amplitude: 

(a\ f [12] (13)^ [46]^ , (12)[26]^(35)^ ] 

M(12;ci_(3)..(4),_(5)..(6)) = { [45] [56]..,e(3(4 + 5)6] + (34) (45).345(3(4 + 5)6] I " 

(178) 

The final result is properly symmetric under complex conjugation and interchange of 
3 ^ 6, 4 ^ 5. Notice that the separate terms of the amplitude are singular not only 
at physical singularities where denominators vanishes, such as S45 = 0, but also on 
the plane where (3(4 + 5)6] = 0. Fortunately, this latter singularity can be shown to 
cancel between the two terms of the expression. 

The two remaining amplitudes can be evaluated by the same technique. For the 
remaining Wuggd amplitude, shifting on 4 and 5 gives 

M(12;ci_(3).-(%.(5)^46)) - / (12)[46][2(5 + 6)4)^ [12][35](1(3^^ _ 
I , K )y \ )y+y ) +y )) \ (45)(56)s456[3(4 + 5)6) [34][45]S345[3(4 + 5)6)J 

(179) 

For the 4-fermion amplitude, shifting on 3 and 6 gives 

I (56)s356[4(5 + 6)3) [56]si23[4(5 + 6)3) J 

(180) 

In this case, one should resist the temptation to shift on the lines 5 and 6. That 
gives a simpler expression, which, however, is not actually correct, because the limit 
2; — )■ 00 in the BCFW construction does not vanish. Alternatively, the 4-fermion 
amplitude can be computed very simply directly from the Feynman diagrams using 
only the technology of Section 2. The result is 



M(12;<(_(3K(4);,_(5K(6)) = _(M&±MM + (31Ml±3)5)M | 

I 5565356 S56S123 



:i8i) 



which can be shown to be equal to (180). 



In principle, one could go further to amplitudes with and any number of 
quarks and gluons. In practice, it is possible to avoid even writing these expres- 
sion explicitly. It is quite straightforward to implement the BCFW recursion in a 
recursive computer algorithm [27] . Then these and higher expressions for amplitudes 
would be generated automatically starting from the MHV and anti-MHV amplitudes. 
Alternatively, these higher point amplitudes are computed in [23)121] . 
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8 Conclusions 

In these lectures, I have given an introduction to a set of tools for computing multi- 
parton tree level amplitudes in QCD. Already at the level of this review, we have seen 
that calculations that would be tremendously difficult by textbook methods become 
accessible, or even easy, using the simplifications of spinor products, color ordering, 
and BCFW recursion. For loop diagrams, even more powerful methods are available, 
as you might glean from the references. I hope that this review will put you on the 
road to a better understanding of QCD that will be useful to you in the era of LHC 
physics. 



A Direct computation of spinor products 

In numerical computations with spinor products, it is best to evaluate the spinor 
products directly from 4-vectors, without first computing vector products. The ad- 
vantage is not only in speed of execution. In the case where two lightlike vectors 
are almost collinear, separated by a small angle 9, the vector product vanishes as 9'^ 
while the spinor product vanishes only as 9^. In such a case, working directly with 
the spinor product avoids round-off error [T5] . 

To evaluate the spinor products of lightlike vectors A and B, first let 



1 ^0 > 

-i ^0 < 

Let 

A+ = A^ + A^ B+ = B^ + B^ . (183) 



{A^-iA^)B+ -{B'-tB')A+ 
(AB) = r]AVB- / , , ; — (184) 



and [AB] = -{{AB))*. 

In this definition, the axis 1 has a preferred role. In compensation, the spinor 
product has a singularity in its phase when one of the vectors A, B approaches the 
— 1 direction. This phase choice cancels out of the squares of matrix elements, so any 
other axis may be used. If the preferred axis is a natural axis of the problem such 
as the beam direction, one will frequently encounter 0/0 in numerical evaluations. If 
the beam axis is the 3 direction, it is useful to choose the 1 direction as the preferred 
axis. Then, working in double precision, excessively small denominators appear only 
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for about 1 point in 10^; in a Monte Carlo integration, one can trap for and ignore 
these points. 
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